MARKOVIAN LOOP CLUSTERS ON GRAPHS 



YVES LE JAN AND SOPHIE LEMAIRE 



Abstract. The loop clusters of Poissonian ensembles of Markov loops on a finite or countable 
^v^j | graph induce a coalescent process on the vertices of the graph. After a description of some 

general properties of the coalescent process, we study the loop clusters defined by a simple 
random walk killed at a constant rate on three different graphs: the integer number line Z, 
the d-dimensional integer lattice Z d with d > 2 and the complete graph. 
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Introduction 



The notion of Poissonian ensembles of Markov loops (loop soups) was introduced by Lawler 
and Werner in [8] in the context of two dimensional Brownian motion (it already appeared 
informally in [13]): the loops of a Brownian loop soup on a domain D € C are the points of a 
Poisson point process with intensity a/i where a is a positive real and \i is the Brownian loop 
measure on D. Loop clusters induced by a Brownian loop soup were used to give a construction 
of the conformal loop ensembles (CLE) in [14] and [12]. They are defined as follows: two loops 
1 £ and £' are said to be in the same cluster if one can find a finite chain of loops £q, . . . , i n such 



that £q = £, £ n = £' and for all i G {1, . . . , n} n £i 7^ 0. Poissonian ensembles of Markov 
loops can also be defined on graphs. They were studied in details in [9] and [6]. In particular 
Poissonian ensembles of loops on the integer lattice Z 2 induced by simple (nearest neighbor) 
random walks give a discrete version of Brownian loop soup (see [7]). 

The aim of this paper is to study loop clusters induced by Markov loop ensembles on graphs. 
The facts presented here, which are built on the results presented in [9] are more elementary 
than in the Brownian loop soup theory but point out that Poissonian ensembles of Markov loops 
and related partitions are of more general interest. The examples we develop show relations 
with several theories: coalescence, percolation and random graphs. 

In section 1, we recall the different objects needed to define the Markovian loop clusters on 
a graph and the coalescent process induced by the partition of vertices they define. In section 
2, we state general properties of the clusters and establish some formulae useful to study the 
semigroup of the coalescent process. In the last three sections we address several aspects of the 
loop clusters induced by a simple random walk killed at a constant rate k on different graphs: 

• On the integer number line Z, the loop cluster model reduces to a renewal process. 
We establish a convergence result by rescaling space by y/e, killing rate by e and let e 
converge to 0. 
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• The loop cluster model can be seen as a percolation model with two parameters a and 
k. Bernoulli percolation appears as a limit if a and k tend to +00 so that converges 
to a positive real. We show that a non-trivial percolation threshold occurs for the loop 
clusters model on the integer lattice Z d with d > 2. 

• As a last example, we consider the complete graph. We give a simple construction 
of the coalescent process on the complete graph and we deduce a construction of a 
coalescent process on the interval [0, 1]. Letting the size of the graph increase to +00, 
we determine the asymptotic distribution of the coalescence time: it appears to be 
essentially the same as in the Erdos-Renyi random graph model (see [4]) though the 
lack of independence makes the proof significantly more difficult. 

1. Setting 

We consider a finite or countable simple graph Q = (X,£) endowed with positive con- 
ductances C e , e G £ and a positive measure k = {k x , x G X} (called the killing mea- 
sure). We denote by C X)V the conductance of the edge between vertices x and y and set 
\ x = ^2 ye x {x y}e£ C x ,y + K x for every x G X. The conductances and the killing measure k 
induce a sub-stochastic matrix P: P x>y = I{{, Ejy } g £} Vx, y G X. We assume that P is 
irreducible. 

1.1. Loop measure. A discrete based loop i of length n G N* on Q is defined as an element 
of X n which can be extended to an infinite periodic sequence. Two based loops of length n 
are said equivalent if they only differ by a shift of their coefficients. A discrete loop is an 
equivalent class of based loops for this equivalent relation (an example is drawn in Figure 1). 
Let VC(X) (resp. VC{X)) denote the set of discrete loops (resp. discrete based loops) of 
length at least 2 on X. We associate to each based loop t = (x±, . . . , x n ) of length n > 2 the 
weight fi(£) = \Pxi,x 2 Px2,x 3 ■ ■ -Px n ,xi- This defines a measure /i on V£(X) which is invariant 
by the shift and therefore induces a measure ^ on VC{X). 





Representation of a loop I of length 4 on the com- 




plete graph K4 corresponding to the equivalence 




class of the based loop I = (4, 1, 2, 1); the measure 




of this loop is 




H(£) = /'l/'ij/'-.i /':.;• 



Figure 1. 



Remark 1.1 (Doob's h-transform) . If h is a positive function on X such that {P — I)h < 0, 
a new set of conductances and a new killing measure k** 1 * can be defined as follows: 

C x h y = h(x)h(y)C Xj y for every {x,y} G £ and n x h ^ = h(x)[(I — P)h](x)X x for every x G X. 
This modification corresponds to the Doob's /i-transform: the associated transition matrix 
verifies P x h y = ^AP x<y Vx, y G X. It is a self-adjoint operator on L 2 (/i 2 A) since = T^PT^ 
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L 2 (/j 2 A) — > L 2 (A) 

with Th : . The loop measure a is invariant under Doob's /i-transform. 

/ ^ hf 

The loop measure can be defined without assuming that P is a substochastic matrix. A weaker 
condition would be to assume that P is a positive and contractive matrix. Taking such matrices 
does not extend the set of loop measures. Indeed, Perron-Frobenius theorem states that if Q is 
a positive and irreductible matrix and if its spectral radius p{Q) is smaller or equal to 1, then 
there exists a positive function h on X such that Ph = p{Q)h. The /i-transform of Q is then 
a substochastic matrix. 

1.2. Poisson loop sets. Let T>V be a Poisson point process with intensity Leb<8>^ defined 
on R + ®VC{X). For a > 0, let VC a denote the projection of the set VP n ([0, a] x VC{X)) 
on VC(X); (VC a ) a >o is an increasing family of loop sets which has stationary independent 
increments. It coincides with the family of non-trivial discrete loop sets induced by the Poisson 
Point process of continuous-time loops defined in [9] . 

1.3. Coalescent process. An edge e G 8 is said to be open at time a if e is traversed by at 
least one loop of VC a . The set of open edges defines a subgraph Q a with vertices X. The 
connected components of Q a define a partition of X denoted by C a . The elements of the 
partition C a are the loop clusters defined by VC a (an example is drawn in Figure 2). This 
paper is devoted to the study of C a . 



C a is a partition with 4 blocks: 
-> {1,2,5,6}, {3,4,7}, 
{8, 10, 11, 12}, {9}. 



FIGURE 2. Example of a loop set T>C a at time a on a finite graph Q, the 
subgraph Q a defined by the open edges and the partition C a induced by its 
loop clusters (the dotted- lines on the left figure represent the edges of Q; VC a 
consists of three loops of length 2, two loops of length 3 and two loops of 
length 4.) 




Remark 1.2. If A is a subset of X, let us define the I?i2 Q -neighborhood of A: 
T a (A) = AU {x e X, 3£ E VC a visiting A and x}. 

Given any (x,y) £ X 2 , set x ~ y if and only if 3k G N* such that y G T a{{ x }) (in the 

a 

example drawn figure 2, r Q ({10}) = {10,11} and for every k > 2, r^({10}) = ^({10,11}) = 
{8, 10, 11, 12} for instance). This defines an equivalence relation and the associated partition 
is C a - 
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2. General properties of discrete loop clusters 

2.1. The distribution of the set of primitive discrete loops. A discrete loop £ is said 
to be primitive if it does not have a non-trivial period; Any discrete loop £ € T>£(X) can be 
represented as some multiple [r/] A/ ^^ of a primitive discrete loop r\ denoted by 7r£. Let VC{X) 
denote the set of primitive discrete loops X of length at least 2 and let VC a denote the set of 
primitive discrete loops defined by VC a . Clearly C a depends only on VC a . 

Proposition 2.1. The distribution ofPC a is a product measure u on {0, 1} VC W defined by 

u(u, uj rt = l) = \-{\-ii{r ] )) a . 

Proof. If r]i,...,r) r are distinct primitive discrete loops, the sets Li = {£ G VC(X), tt£ = rji}, 
i = 1, . . . , r are disjoint, hence the r events Ei = {3£ 6 T>£ a such that ir£ = r/i}, i = 1, . . . , r 
are independent. Therefore, the law of VC a is a product measure on {0, 1} VC W and for every 

+oo 

P(r] e VC a ) = F(3£ G tt! = r?) = 1 - exp(-a ^ fi(£,t= [vD)- 

m=l 

As Ai([??] m ) = ^fi(r]) m for every m € N*, we deduce that 

P(t? e P£ Q ) = 1 - exp(alog(l - (i(rj))). 

□ 

It follows from Proposition 2.1 that Harris inequality (and also B-K inequality) holds on 
increasing events (see e.g. [15]). In particular let us say that a subset A of X is connected at 
time a if it is contained in a cluster of C a . Then by Harris inequality, 

¥(A and B are both connected at time a) 

> ¥(A is connected at time a) ¥(B is connected at time a). 

In a similar way, we say that a subgraph is open at time a if all its edges are traversed by jumps 
of loops in T>C a . Then the same inequality holds. 

These inequalities can be extended to any number of increasing events. In particular, 

• For every partition ir = (Bj)j 6 / of X into non-empty blocks, the probability that C a is 
a coarser partition than ir satisfies: 

P(C Q H 7r) > JJP(-Bj is connected at time a). 

i&I 

• If J- C £, then P(J 7 is open at time a) > TT P(e is open at time a). 
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2.2. The transition rate of the coalescent process. The evolution in a of C a defines a 
partition- valued Markov chain. Let tt be a partition of X into non-empty blocks {Bj, i G I}. 
From state tt, the only possible transitions are to a partition 7r®^ obtained by merging blocks 
indexed by some subset J of / to form one block Bj = Uj^jBj and leaving all other blocks 
unchanged. The transition rate from tt to tt® j is 

r n n ®j = [/,(£, Vj £ J, £ intersects Bj and Vu £ J, £ does not intersect 2? u ) 

= T ^! -fxi,X2 ■ ■ ■ -fajfc,a:i (2-1) 

fc>2 (xi,...,a! fc )ew fc (^,ieJ) 

where Wk{Bj, j £ J) is the set of £>tuples of elements of Bj having at least one element that 
belongs to each block Bj, j G J: 

Wk(Bj, j £ J) = {(x\, . . . ,Xk) G (Uj£jBj) k such that Vj £ J, x u £ Bj for some it}. 

If the graph Q is finite, these transition rates can be expressed with determinants of Green 
functions of subgraphs. In order to describe the formula, let us introduce some notations; Let 
G denote the Green function of Q: G = (XI — C)" 1 where A = (X x ) x ^x-, I is the identity matrix 
and C is the conductance matrix. For every subset F of vertices, set F c = X \F and let G^ 
denote the Green function of the subgraph (F,£\ FxF ) of Q: G^ F > is the inverse of the matrix 
(XI - C)\ FxF . 

Proposition 2.2. Let us assume that Q = (X,£) is a finite graph. 

For every partition tt = {Bj, i £ 1} of X and every subset J of I with at least two elements: 

T v>n *J = Et" 1 )'' 1 log(det(G( u «"\A))). (2.2) 

Proof. By the inclusion-exclusion principle, 

t,. _ej = u[£, £ does not intersect U Bj) — till, 3u £ J, I does not intersect U Bj) 
V j£J V V ieJ c u{n} v 



E (~ l ) lK ^(^ t does not intersect B u ). 



To conclude, we use that: 

• For a subset F of vertices, 

H(£, £ does not intersect F) = fi(VC(F c )) = log ( det(G (i?c) ) ] [ X 3 

x£F c 

• For every family of reals (a u ) ue j indexed by a finite subset J with at least two elements, 

Em) 1 " E a « = °- ( 2 - 3 ) 

/cj ueJ\i 

□ 

Example 2.3. Let us consider the complete graph K n with unit conductances and a uniform 
killing measure with intensity k. The transition matrix P verifies: P XjV = n _\ +K f° r ev_ 

ery vertices x and y. Thus for every subset F of vertices, det(G^) = ((n + K)l i? l _1 (n + k — ji 7 ))) -1 
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(the computation of this determinant is detailed in Lemma A. 3). Using (2.2) and (2.3), we 
obtain: 

v^Ec-^+M 1 -^ £ M- (2 - 4) 

/gj ueJ\-f 

2.3. The semigroup of the coalescent process on a finite graph. In this section, we 
assume that the graph Q is finite. The probability that C a is finer than a given partition of X 
has a simple expression: 

Lemma 2.4. Let us assume that the graph Q = (X, £) is finite. Given two partitions ttq and 
TT = (Bi) ieI of X, 

r„ [C a b ir) - I det(G) J {7rob7r} ' ^ ' 

Proof. Let us assume that tt is coarser than ttq. The event 'C a is finer than tt' means that 
every loop of I>£a is included in a block of tt. Therefore, 

K (C a h tt) = exp ( - a(ji(VC(X)) - ^//(P£(^))^ 

since the set of loops in each block i?j of 7r at time a defines independent Poisson point processes 
with intensity pS- Bi ^ respectively. To conclude, we use that 

H{VC{X)) = log(det(G) X x ) and n{VC(Bi)) = log(det(G^) X x ). 

□ 

An explicit formula for ¥ no (C a = tt) can be derived from Lemma 2.4. Let us first introduce 
some notations. For a partition tt, let \tt\ denote the number of non-empty blocks of tt. For a 
subset A, let tt\a denote the restriction of tt to A: tt\^ is a partition of A, the blocks of which 
are the intersection of the blocks of tt with A. 

Lemma 2.5. Let us consider a finite graph Q = (X,£). Let ttq and tt be two partitions of X. 
If tt has k non-empty blocks denoted by B\ . . . ,B^ then 

k 

^ (C a = tt) = EH^IIO^J " !) !IP -o(C a h tt). (2.6) 

7r>^7r i=l 

Proof. Let us first assume that tt = {X}. To obtain equation (2.6), it is sufficient to prove the 
following identity: 

\x\ 

%.=w} = E(- 1 )'" 1 (^- 1 ) ! E %^>- ( 2 - 7 ) 

where ^3g(X) denotes the set of partitions of X with £ non-empty blocks. Let us assume that 
C a is a partition with j non-empty blocks. For I < j, we can construct a partition coarser than 
C a with I non-empty blocks by choosing how to merge some blocks of C a , that is by choosing 
a partition of {1, ... , j} with I blocks. Therefore 

E 1 {c a > n } = m({i,---j})\i { i< j} 
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and the right-hand side of (2.7) is equal to J2e=i(~ 1)^ _1 (^ — 1)'|^({1> • • • > By an identity 
on the Stirling numbers of the second kind (see for example [11] equation (1.30) page 22), this 
sum is equal to 1 if j = 1 and if j > 2. 

To prove (2.6) for a partition 7r with A: non-empty blocks {B\, . . . , B^), we consider C Q as a 
sequence of partitions on Bi, . . . , respectively and apply (2.7): 

k \Bi\ \B k \ k 

I{C=(Bi,...,B fc )} = n%«| Si =B»} = E " EII^' 
i=l £i = l 4=1 i=l 

where = (-F/" 1 ^- l)!E#e%(B0 ^C,,^.^}- Therefore, 
|Bi| \B k \ 

t { c a = {Bl ,...,B k)} = E-- E(- 1 )' 1+ - +4 ""(^- 1 ) ! ---(4-i)! E 

4=1 4=1 7fe% 1 (Bi)x...x% fc (B fc ) 

7r^7r i=l 

□ 

Example 2.6. Let us consider the complete graph K n endowed with unit conductances and a 
uniform killing measure of intensity k. If tt is a partition of the set of vertices X with k blocks 
B\ , . . . , Bk then 

p,„(c„ t ,) = (^rn(i-JfL)-°W.)- 



and 



P. (C Q = n) = (_) E (-l)' Jhfc n ( card (^ G J ' ^ C n (i- — 

S.t. 7TQ >^ 7T >^ 7T 

Let us note that (1 7— ) _a is the 7-th moment mi of the random variable 1" = expf^f— ) where 

Z denotes a Gamma(a, l)-distributed random variable. Thus P no (C a >z it) = — ^ I{7r ^7r}- 

Ut=i m \Bi\ 

Let c n denote the n-th cumulant of Y. Formula (2.6) and the expression of cumulants in terms 
of moments (see formula (1.30) in [11] for example) yield ¥ no (C a = {X}) = 

2.4. Loop clusters included in a subset of X. For a subset D of X, let VC^ denote the 
loops of T>C a contained in D and the partition of D induced by T>C^ . In general 
is finer than the restriction of C a to D but coincides with it if D is a union of elements of C a . 
For every partition it of X with k non-empty blocks Bi, . . . , B^, the loop sets VC^ 1 ^, . . . , 

(B ) h (B ) 

VC y a , VC a \ (Llf =l VCa ) are independent. This yields the following equality: 

k 

F(C a = vr) = F(C a £ tt) JJP(4 Bi) = {Bi}). 

i=l 

Let us note also that if U CflCl then 

>- {U, U c }) = F(Ci D) h{U,D\ U}) F($£ e VC a visiting U and D c ) 
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(see [5] page 2 for a related result in Schramm-Loewner Evolution (SLE) context). The second 
term on the right-hand side equals 

e -a(v(VC(X)y^VC(U c ))-v(VC(D)) +lJ ,(VC(D\U))). 

If X is finite, it can be written as 

det(G( l/c ))det(G( D ))x" _ ,., f det(G ]UxU )det(G ]DCxD o)^-a 

or, with Jacobi s identity, as 1 



det(G)det(G( D W)) J ' " V det(G| (C7uD c )x(C7uDC) ) 

2.5. Computation of the semigroup using exit distributions. The probability that C a is 
finer than a partition tt can also be expressed using exit distributions. The formula obtained is 
easier to use than (2.5) for graphs such as tori, trees . . . Let us first introduce some notations. 
For a subset D of X, let dD denote the inner boundary of D: 

dD = {x G D, C X)V > for some y G D c } 

and let denote the exit distribution (or Poisson kernel) from D: for x G X and y G D c , 

H^j is the probability that a Markov chain with transition matrix P starting from x exits 
from D at y. 

Lemma 2.7. Let tt = (-61)1=1,...,*; be a partition of X. Let B = Uj =1 dBj denote the union of 
the boundary points of the blocks of tt. Let %M denote the matrix indexed by B defined by: 



1 {x= y } ifx,yedBi 
-Hxfy^ if x £ dBi,y G dBj and i ^ j. 

IfB is finite then F(C a £ tt) = det(^W) Q . 



Proof. 



• First, let us assume that X is finite. Let K denote the product of the block diagonal 
matrix diag(G^), i G {1, . . . , k}) by the matrix G~ l . We can rewrite the expression of 
F(C a £ tt) given by Lemma 2.4 as P(C a £ vr) = det(K) a . The restriction of K to B { x B t 
is the identity. The exit distribution from a subset D verifies: Hj^y = ^2 z ^d G^)c z ^ y 
for every x G D and y G D c . Therefore, 

KC t ») - *W with if., = { *< 5 |, ^ £ ^ md ( # . (2.8) 

Let (£ n )n denote a Markov chain with transition matrix P. The trace of (£ n )n on B 
defines a Markov chain denoted by (£ n ) on £> and thus a Poisson point process T>V on 
R + xVC(B) (see [9] chap. 7): the discrete loops set at time a is 

VC a :={l\ B , leVC a \v6* c) }- 

Let C a be the partition of B induced by VC a . The non-empty subsets dBi, i G {1, . . . , k} 
define a partition of B denoted by dir. As {C a £ it} is the event 'Vi, j G {1,. . . ,k} 
such that i 7^ j, Vx G dBi, y G dBj, {x,y} is not crossed by a loop ofVC a at time 
a', it only depends on the restriction of the loops on B, hence P(C a £ tt) = P(C a £ dir). 
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For a subset F of B, let H^ F > denote the exit distribution from F for (£ n ). It follows 
from formula (2.8) applied to C a that 



F(C a h tt) = det(K) a where K x , 



1 {x= y } i f x, y G dB,, 
H x d y l) if x G dB u y G ftB, and i / j. 



~ (QB- ) ( B ■ ) 

To conclude, it remains to note that H XtV = H Xj y for every x G dBi and y G B\dB{. 
• Let us now assume that X is a countable set. Let (X^^ be an increasing sequence 

of finite subsets of X such that X = U^^Xk- As B is assumed to be finite, there 

( x 

exists an integer &o such that B is included in Xk (r For k > kg, a loop in "D£a that 
passes through two different blocks of tt, passes through two different blocks of TT\x k 
the restriction of tt to X^. Therefore, 

F(C a >z tt) = 1 - sup F(3£ G VC^ k) passing through two different blocks of tt) 

k>ko 

= wP(C«>t^). 

By the first part of the proof, P(ci Xfc) >: TT\ Xk ) = det(H {n * x k ) ) a . It remains to 
note that the matrix T-L^ x k' coincides with for every k > ko to deduce that 

h tt) = det(H^) a . 



□ 

2.6. Closed edges in a finite graph. An edge is said to be closed at time a if it is not 
crossed by a loop of VC a . More generally, a family of edges E is said to be closed at time a if 
every edge of E is closed at time a. 

Let us assume that Q is a finite graph. To compute the probability for a family of edges E 
to be closed, we modify the conductances and the killing measure so that the conductance 
of every edge in E is and the measure A is unchanged: C XjV = C XjV ^{{ x , y }^E} G £ 

and k x = k x + J2 y { x y}eE &x,y Vx G X. Let Ge denote the Green function associated with 
{C e ,e G £, k x ,x G X}. 

Lemma 2.8. Let us assume that Q = (X,£) is a finite graph. The probability for a family of 
edges E to be closed at time a is ■ 

Proof. For a loop £ and an edge e, let N e (£) denote the number of jumps of £ across e. Set 

P(E is closed at time a) = exp ( - atfi(VC(X)) + afi(£, N E (£) = 0)) . 

Recall that fi(VC(X)) = — log(det(I — P)). For s G [0,1], let Pe(s) denote a perturbation of 
P defined by 

\ sP xw if {x, y} G E 
\ P x , y if {x,y} $E. 
Similarly, [i(s Ne ) = — log(det(I — Pe(s))). In particular, 



[Pe(s)] xa 



H(£, N E {£) = 0) = -log(det(J - P E (0))). 
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Lemma 2.8 follows since Pe(0) is the transition matrix associated with the {C e , e S £, k x ,x £ X}. 

□ 

We can deduce from Lemmas 2.8 and 2.4 that if E is a set of edges of a finite graph with 
extremities in different blocks of a partition tt = {L>i, i £ 1} then 



P(C Q y 7T I E is closed at time a) = ( 



n 467 det(G<(^ 



det(G s ) y 

3. Renewal processes 

On the graph Z, the clusters C Q are intervals between closed edges at time a (namely edges 
which are not crossed by any loop of T>C a ). The Poisson loop sets induced by a simple random 
walk killed at constant rate k have the following properties: 

Proposition 3.1. Let us consider the graph Z endowed with unit conductances and a uniform 
killing measure with intensity k. Set = log(l + § + \J k + ^-). 

(i) The midpoints of the closed edges at time a form a renewal process. Moreover, 

— the probability that {n,n + l} is closed at time a is equal to (1 — e _2pW ) Q for every 
n e Z; 

— given that {0, 1} is closed at time a, the probability that {n, n + 1} is also closed 
is equal to 1 (k)( — J for every n G N. 

(ii) Assume that a g]0, 1[. Let denote the law of this renewal process i.e. the law of 
the distance between two consecutive closed edges at time a. For e > let (^ £ ,i)ieN* 
denote a sequence of independent random variables with distribution i/( K£ ) . 

l-a 

[e~ J-t] 

For every t > 0, as e converges to 0, y/e Y e ^ converges in law to the value at t 



of a subordinator with potential density {- 



i=l 

2y^ 



Remark 3.2. This subordinator is associated with the Poisson covering (cf. [1], chapter 7) 
defined by the infinite measure on M + with density u i— > ~ a ^p: log(l — e~ 2u ^) = ^ K " ( L 2 Tj^2 ■ 
These covering intervals can be viewed as images of the "loop soup" of intensity a associated 
with the Brownain motion killed at constant rate k. 

Proof of Proposition 3.1. 

(i) An edge {x — l,x} is closed at time a if and only if VC a = VCa +N ^ UVCa 1 . 
The next closed edge is the first edge which is not crossed by any loop of VC a , 
and previous closed edges are defined by VCa 1 N ^ which is independent of VCa +N ^ . 
Stationarity is obvious as VC^^ — x is distributed like VC^. For n G N, let 
r^ denote the probability that {n, n + 1} is closed at time a. By Lemma 2.7, 
ri^ = (1 — H^^^^H^^ ) a where Hx°y denote the probability that a simple ran- 
dom walk killed at constant rate k and starting at x exits D at point y. The functions 
x i — y can be computed from the solutions of equation 

(2 + n)u{x) - u(x + 1) - u(x - 1) = 0, 
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which are linear combinations of exp(p( K ) x) and exp(— p^x). 

We obtain H%+£ N) = H^~^ kjn+1 = exp(- k) for every k G N. 

For n G N*, let q( K \n) denote the probability that {n, n + 1} is closed at time a given 

that {0, 1} is closed at time a. To compute q^ K \n), we consider a simple random walk 

on N killed at rate k and at point 0; we denote it (Cfc)fc- As (V£^) a is the Poisson loop 

sets associated with (Ck)k, by Lemma 2.7 we obtain: q^ K \n) = (1 — /i^ °''"' n ^/i^ 1 1+N ^) 

where hS D > is the exit distribution from a subset D of N for (Ck)k- 

= Sj ; ih( f Wfc) ^ for 1< fc < n + 1 and 4"+ 1+N) = for A; > n. 

fc sinh(p(»)(ra + 1)) - - fe 



Therefore a^(n) = (1 sinh( P («)n) - p (*K a = / i- e - 2 P W 

inereioreg in; ^ sinh(p («) (n+ i)) e J S_ e -2 P («)(n+i) 



(ii) Let I re denote the Laplace transform of the function u i— )■ ( ^y^ )"- 

/ K ( s ) = ^ O °( T -^re-^for S >0. 

Let i^ K ) denote the Laplace transform of z/ K ). We shall prove that for every s > 0, 
— e^2~ log(z/ K£ ) (si^/ej) converges towards l/I K (s) as e tends to which yields (ii). 
As q( K \n) = ^2'f? = i{v^)* k (n) for every n £ N*, we have the Laplace transforms iden- 
tity: z/( K ) = L . Therefore, it is sufficient to show that e^ 1 ~ Q ^ 2 g , ( K£ )(sy / e) converges 

1+qM 

to 7„(s) as e tends to 0. Let n\ and «2 two positive reals such that < K\ < k < k%. 
For e small enough, sJek\ < p^ K£ ^ < \j£K2 and 2y/£H\ < 1 — e _2p ^ < I^Jeki- For 
a > let / a be the function defined by f a {x) = (1 — e _2aa; ) _a for x > 0. Therefore, for 
e small enough, 

(2y^T) a / K2 ( v / F(7i + l)) <q (Ke) (n) < (2^) a f Kl (yfe~(n + l)) for every n G N . 

By Lebesgue's dominated convergence theorem applied to the function g a , s ,e defined by 
9a,s,e{x) = f a (^\^\)e~ sV ~ £[ ^ Vx > 0, we obtain: 

V^y2fa(V^(n + l))e- s ^ n -> / f a {x)e- sx dx 

(for e < i, 5 a , Sj£ is dominated by i ^ (^)° 1 ]0i _lj + {1 _ e e -7/2 )a I]-L ;+oo[ (^))- This 
implies that for every s > 

+°° ' 



o 



1 _ g-2 v ^in 



su du < lim e^-^q^sy/e) 

<BEe( 1 - ot )/ 2 ^)( s Vi) < /^f Ye~ su du 

Jo VI - e" 2 ^"/ 



These inequalities hold for every < ki < k < «2- Therefore for every s > 0, 
£( 1-a )/ 2 g(Ke) converges to I K (s) which ends the proof of (ii). 



□ 
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Remark 3.3. In the case of the simple random walk on N killed at 0, a similar argument 
(detailed in [10]) shows that: 

• For < a < 1, the midpoints of the closed edges at time a form a renewal process 
with holding times (Yn)n>i\ The generating function of is 1 — Li s ^ where Li 

denotes the polylogarithm: V |s| < 1, U a (s) = Y$=i Set S ^ = E?=i for 
n > 1. As e tends to 0, { £ S^l-i t \ > t ^ 0) converges in law towards a stable subordinator 

(sj: a \t > 0) with index 1 — a. In the case of a finite interval [0, L], we obtain a renewal 
process conditioned to jump at point L. 

• For a > 1, there are only a finite number of clusters. In particular, P(S'| Q ^ = oo) = 

4. Bernoulli percolation and loop percolation 

4.1. Bernoulli percolation. Let s = {s e , e £ £} be a family of coefficients in [0,1]. In 
Bernoulli percolation model of parameter s on the graph Q = (X,£), every edge e is, inde- 
pendently of each other, called 'open' with probability s e and 'closed' with probability 1 — s e . 
Vertices connected by open paths define a partition of X denoted by V(s). We can compare 
V(s) to the partition induced by the set of primitive discrete loops of length 2 for a finite graph 
(X,£). 

Lemma 4.1. Let us consider a finite graph (X,£). For a > and e = {x,y} £ £, set 

Sa,e — 1 (1 Px,yPy,x) ■ 

The partition induced by the set of primitive discrete loops of length 2 at time a has the same 
law as V({s ate , e € £}). 

Proof. Let VA a denote the set of primitive discrete loops of length 2 at time a. The law VA a 
for every a is a product measure v on {0,1} £ : for every edge {x,y} S £, let rj Xi y denote the 
class of the based loop (x,y). The probability that a loop in VC a has {x,y} as support is 

u(oj, LJ{ x , y } = 1) = 1 - exp ( - afi(£, ir£ = r) XjV )) . 

and 

+oo ^ 

M(A = Vx,y) = Yl ^(Px,yPy,x) k = -log(l - P x ,yP y , x ). 
k=l 

□ 

Bernoulli percolation clusters on a graph appear to be a limiting case of partitions defined 
by loop clusters in which only two points loops contribute asymptotically. 

Proposition 4.2 ([10]). Let us consider a finite graph Q = (X,£) endowed with unit conduc- 
ts 

tances and a uniform killing measure with intensity k > 0. Let Ca be the partition induced by 
the Poisson loop set on Q at time a. Fix u > 0. 

If k and a tend to +oo such that a is equivalent to uk 2 , then converges in law towards the 
Bernoulli percolation of parameter 1 — e~ u . 
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Proof. For a partition it = {Bi, i G 1} of X, let F(7r) denote the set of edges of £ linking 
different blocks of ir. 

The law of Bernoulli percolation of parameter 1 — e~ u denoted by V(l — e~ u ) is characterized 
by the identities: 

F(V(l - e~ u ) y tt) = e ~ u|i(7r)l for every partition vr of X. 

To prove the convergence of F(C^ y n), we apply Lemma 2.7: P(ci y tt) = det(H^) a 
where is defined as in Lemma 2.7. We note then that Tix^y is equivalent to k~ 1 if {x,y} 

belongs to L(ir) and of order less or equal to k~ 2 otherwise. Indeed, if x G dBi and y G <9-Bj 
for i / j, = £ z£B . P*,*^ 5 + P^ and P hJ < I for all (t,i) eXx! 

A second-order Taylor expansion shows that log(det('H^ 7r ^)) = Tr (log is equivalent to 

-\Tr{Q 2 ), With = At -1 I{{a;,y}eL(7r)}- D 

4.2. Loop percolation on Z rf with d > 2. Let us consider the Poisson loop process induced 
by the simple random walk on Z d , d > 2, killed at a constant rate k > 0: P x , x +u = 23Tk ^ or 
every x G Z rf and u G {±l} rf . Let 6(ol,k) denote the probability of percolation at time a i.e. 
the probability of any fixed point to be connected to infinity by an open path at time a. The 
following Proposition presents some properties of the function (a, k) i— >■ 9(a,n): 

Proposition 4.3. Let p c denote the critical probability for bond percolation on Z rf (d > 2). 

(i) #(a,K;) is an increasing function of a and a decreasing function of k. 

(ii) 9(a, k) > for every a > and k > suc/i that (1 — ( 2 d+ K )^ )" < 1 — -Pc- 

(iii) For any a > 0, 6(a, k) vanishes for n large enough. 

Proof. 

(i) 9(a, k) is an increasing function of a since a i— )■ T>Ca is increasing. 

To show that «) is a decreasing function of k, we use an independent thinning 
procedure. Let k\ > K2 > 0. The corresponding measures on based loops satisfy: 
= = i(2d^) fcl {x I+1 -x ie {±i^Vie{i,...,fc}} f or j G {1,2}, fc > 2 and 

xi, . . . , Xk G Z d with the convention x^+i = sci. By erasing independently each based 
loop I G VC^a 2 ^ of length > 2 with probability 1 — ( f§j^ ) fc , we obtain a discrete loop 
set having the same distribution as VC^ 1 ^. 

(ii) By Lemma 4.1, the partition induced by the set of primitive discrete loops of length 2 at 
time a has the same law as the Bernoulli percolation with parameter 1 — (1 — ^d+n)' 1 

It follows from Bernoulli percolation on Z d that if 1 — (1 — ^d+n)' 2 )° ^ ^ c then 
6>(a, k) > 0. 

(iii) To prove that 9(a, k) vanishes for k large enough, we show that there exists a finite 
real C a > such that any self-avoiding path x = (x±,X2, ...,x£) of length L G 2N* 
is open at time a with probability less than (%0 L - We can then conclude with the 
usual path-counting argument: for every L G 2N*, 9(a,n) is bounded above by the 
probability that there exists an open self-avoiding path of length L at time a starting 
from the origin, hence 9(a,n) < limsup L _ s>+00 (2d) L (^) L = for k > 2dC' a . 

Let us first introduce some notations. 
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— Let ^3(2, {1, . . . , L}) consist of partitions of {1, . . . , L} in which all blocks have at 
least two elements (the number of blocks of such a partition ir is denoted by |7r| 
and blocks are denoted by 7Ti, ir 2 , ...). 

— For a vertex v and a loop £ G T>£(Z d ), let N v {£) denote the number of times £ 
passes through v. N v (£) = \{1 < i < k, Ui = v} | if £ is the class of the based loop 

(m, . . . ,u k ). 

Let us consider a self-avoiding path of length L G 2N* denoted by x = (x\,X2, . . . ,xl) 
and let £ x denote the set of edges {x 2 i-i,x 2 i}-, 1 < i < L/2. 

To be open, the edges of x have to be covered by the edges of N < L loops of VC a . 
Among these loops, those that cover at least one edge e € £ x can be used to define a 
partition ir G ^3(2, {1, . . . , L}) as follows: let £\ be a loop in VC a covering edge {x±, x-i\. 
The first block ix\ of ir consists of the indices of the endpoints of e G £ x covered by £\. 
If 7Tx 7^ {1, . . . , L}, let j be the smallest integer i such that xi is not in tti and let £2 
be a loop covering {x^, The second block ir 2 of 7r is defined as the set of indices 

of endpoints of e G £ x covered by £2 that are not in tti, and so on until all elements 
of {1, . . . , L} belong to a block written down (an example is presented in Figure 3). 
Therefore, 

M 

I{z open at time a} < ( II II (^)>°> ) 

tt6 ( P(2,{1,...,L}) <?i,...,<?| 7r |eX>£ Q j=li€irj 
pairwisc distinct 

s e ( e nn^ft))- 

7re s P(2,{l,...,L}) £i,.../ K |eX)£ Q j=l tSTTj 
pairwisc distinct 



Example of an open path x = (1,2,3,4,5,6,7,8) 
with five loops covering its edges. Three of them 
(drawn in solid lines) cover the set of edges £ x \ 
we associate to these three loops, the partition 
vr G 93(2, {1, . . . , 8}) with three blocks = {1, 2}, 
vr 2 = {3,4,7,8} and ^3 = {5,6}. 

Figure 3. 

By Campbell formula (equation (4.2) page 36 in [9]), for every k G N* and for every 
positive functions F±,...,Fk 

e ( e nw) = n ( / 

h,...,l k eDC a , i=l i=l 
pairwisc distinct 

Therefore, 

kl , 

P(x is open at time a) < ^ [ / ( ] I N Xi {£)^a^{d£). 
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For a finite set A, let &a denote the set of permutations of elements of A. By Propo- 
sition 5 page 20 in [9], for every k > 2 and for every vertices yi, ■ ■ ■ ,yk, 

i=1 i=1 ^6{i,...,k} 

L 

Therefore, P(x is open at time a) < ^ JJ A^J S a (x) where 

i=l 

7re«P(2,{l,...,L}) j=l 1 ■ ?l <jG6^ 

The blocks of a partition -zr 6 ^3(2, {1, . . . , L}) in (4.1) can be seen as the orbits of a 
permutation without fixed point. Since a circular order on k integers corresponds to k 
different permutations of these integers, S a (x) can be rewritten as follows: 

S a (x)= £ a m(CT) G^ CT(1) ...G, iW (4.2) 

where 

— m(a) denotes the number of cycles in a permutation a, 

— L i denotes the set of permutations of {1, . . . , L} without fixed point (such a 
permutation corresponds to a configuration without isolated points). 

(See Lemma A.l in the appendix for details.) The right-hand side of equality (4.2) is 
nothing but Per° (G XijXj , 1 < i,j < L) defined in [9] page 41. 
To conclude we use that 

L 

Per?(G w l<i,j<L)<;Q( £ G Xi , Xj ) < ( £ Gq,,) 1 ' 

»=1 l<j<L, j^i yf zZ d \{Q} 

since the vertices xi are pairwise distinct. As (Pl) u = for every u G Z , 

£ A,G 0>y = £(P fe l)o = ^. 
s/ez d \{o} fc=i 

Thus P(rr is open at time a) < (^ max(a, 1)) L which ends the proof. 



□ 



Remark 4.4. 



(i) It follows from Proposition 4.3 that for every a > 0, there is a finite value of k, we can 
denote n c (a), above which 6(a, k) vanishes; Morever k c is an increasing function that 
converges to +oo as a — > oo. 

(ii) As the simple random walk in Z 2 is recurrent, for d = 2 the probability that a fixed 
edge is open at time a > converges to 1 as k tends to 0. Indeed, let us first note that 
for every x G Z 2 , F(N^ a) = 0) = ( Xx g xx Y ( this equality has been proven for a finite 
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graph 1 , so we apply it to the restriction of the random walk to Z 2 n[— M, M] 2 and let 
M tend to +00). Therefore F(N^ a) > 0) converges to 1 as k tends to 0. 
Fix x G Z 2 and u E {±1} 2 . By symmetry, 

P(JVW>0)= P (Ni% v > 0) = 4P(A$ fu > 0). 

t>e{±i} 2 

Since a h-» VC a has independent stationary increments, ^{N^ +u = 0) = ^(N^+l = 0) 
for every n > 0. Therefore, 

P(ivg +M > 0) = 1 - P(A^> = 0)" > 1 - (1 - \ P(iVi a/n) > 0))". (4.3) 

For a fixed e > 0, let us choose n such that (1 — g) n < s. There exists k £ such that for 
k < k £ , P(ivi a/n) = 0) < |. It follows from (4.3) that P(A r iJ +u > 0) > 1 - e. 



5. Complete graph 

This section is devoted to the study of Poisson loop sets on the complete graph K n endowed 
with unit conductances and a uniform killing measure. The set of vertices is identified with 
{1, . . . , n} and the set of partitions of {1, ... , n} is denoted by ^3({1, . . . , n}). The intensity of 
the killing measure is denoted by K n = ne n with e n > 0, hence the coefficients of the transition 
matrix P are: P x>y = ^ffi with A^ n ^ = n — 1 + ne n for every x,y E {1, . . . ,n}. As n will 
vary, the loop set and the partition defined by the loop clusters at time a will be denoted by 
T>C^ and respectively. In a first part, n is fixed; We present another construction of the 
coalesecent process (Ca) a >o and use this construction to define a similar coalescent process 
on the interval [0, 1]. In a second part, we let n tend to +00 and describe the distribution of 
the first time when has no block of size one (cover time) and the first time when (C^) a 
has only one block (coalescent time). 



5.1. Another construction of the coalescent process. 

Proposition 5.1. From state it = {Bi, i G 1} G *P({1, . . . ,n}) ; the only possible transitions 
of {C^) a >o are to partitions tt® j where J = {j±, . . . is a subset of I with L > 2 elements. 
Its transition rate from tt to ir® J is equal to: 

1 k 

k>L v n 1 (h,...,i k )eW k (J)u=i 

= E kn k {£n + l)k E L k J II ( 5 - 2 ) 

fcl+...+fc i =fc 

where Wi-(J) is the set of k -tuples of J in which each element of J appears. 



W a finite graph = (X, £), ¥{N ( x a) = 0) = exp ^-a^(VC{X))-fi(VC(X\{x}))^ = ( 
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(ri) 

Proof of Proposition 5.1. Let us recall the expression of t @J obtained in Example 2.3: 



By expanding the logarithm, we obtain: 

T l,l® J = E LM(n) T iU X]^ -1 )' 7 ' ( E I 5 ' 



To establish 



k 



/,/cj ueJ\-f (ii,..4)eff t (J)n=i 



E M = E H^i ( 5 - 3 ) 



which ends the proof of (5.1), we make use of the exclusion-inclusion principle. Let us con- 
sider the following random experiment: 'we choose uniformly at random k points in the set 
B := Uj^jBj 7 and for j E J, let Ej denote the event 'at least one of the k points falls 

into Bj '. By the inclusion-exclusion principle, the coefficient X^/cj(— l)' 7 ' ^ JZ u eJ\i l^l) di- 
vided by \B\ k is equal to P(n je j£'j). The probability ¥(C\j^jEj) can also be decomposed 
by introducing the events A^j 'the i-th point falls into Bj' for every i £ {l,...,fc} and 
j G J- F(r\j£jEj) = Z)(i lj ...,i fe )6iv fc (J r ) P ( n «=i^«>in)- Tnis corresponds to the right-hand side 
of (5.3) divided by |-B| fc . Equation (5.2) is obtained by rearranging the terms of the product 

uLi\B iu \. □ 

Formula (5.1) of the transition rate has a simple interpretation: if we choose R points uni- 
formly at random in {1, ... ,n} where R has the logarithmic series distribution 2 with parameter 
- the probability that at least one point falls into each block {Bi, i € J} and none out- 
side of Uj & jBj is equal to (— log(l — - 1 | _ 1 ))~ 1 r^ ffiJ . From this remark, we derive a simpler 
construction of (C^) a >o- 

Proposition 5.2. Let us define a ^3({1, • • • ,n})-valued sequence {Ykjk iteratively as follows: 

• Yq is a random variable with values in *P({1, . . . , n}); 

• Let k £ N. Given Y^, 

— we choose an integer R>2 independent of Yq, . . . , Y n , following the distribution 

v = ^k>2 p enk (L+l)* 6k where = - - iAi) - i^T 

— we choose R points U\, . . . , Ur uniformly at random on {1, . . . , n} and indepen- 
dently ofY , ...,Y k ,R. 

— Yfc+i is obtained from by merging blocks of that contain at least one of the 
R points Ui, . . . , Ur. 

Let (Zt)t>o be a Poisson process with intensity (3 £n and independent of {Yk)k- 

The process (Jl[ n ^)t>o defined by IIj = Yz t for every t > is a continuous Markov chain: 



2 The logarithmic scries distribution with parameter < a < 1 is defined as the measure 



Un ~ -log(l-a) 2^fc = l k ° k - 
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The only possible transitions from state tt = (_Bj)j 6 / are to partitions tt® j where J is 
a subset of I with at least two elements and its transition rate from tt to tt® j is r m , 
defined in Proposition 5.1. 

If tt is a partition of {1, . . . , n} with k non-empty blocks {B\, . . . , B^) then 

k 

It 

i=l 



% ,nH^) = (^na-JfL_r^ 



5.2. A similar coalescent process on the interval [0, 1]. The algorithm described in Propo- 
sition 5.2 can be adapted to define a coalescent process (IT^oo on the interval [0, 1] such that 
for every partition tt of [0,1] and every t > 0, P(IL^ >z 7r^) converges to P(Iit ^ 7r) if e n 
converges to e > and the partition tt^ of {1, . . . , n} converges to tt as n tends to +oo. 

Proposition 5.3. For k £ N*, let ^([0, 1]) denote the set of partitions on [0,1] with k 
blocks having a positive Lebesgue measure. Set *p([0, 1]) = U/ c6 ^*^fc([0, 1]). For a partition 
tt = (b{)i£i £ *P([0, 1]) and a subset J of I, let tt® j denote the partition obtained from tt by 
merging blocks bi for i £ J. Let e be a positive real. Let us define a *P([0, l])-valued sequence 
(Y n )n iteratively as follows: 

• Yq is a random variable with values in *P([0, 1]); 

• Let n £ N. Given Y n , 

— we choose an integer R > 2 independent of Yq, . . . , Y n following the distribution 

v = Ek>2 fefcfa^ fc where & = - Mi - iii) - iii- 

— we choose R points U±, . . . , Ur uniformly at random in the interval [0, 1] and in- 
dependently of Yq, . . . ,Y n ,R. 

— Y n+ \ is obtained from Y n by merging blocks of Y n that contain at least one of the 
R points Ui, . . . , Ur. 

Let (Zt)t>o be a Poisson process with intensity (3 e and independent of (Y n ) n . 

The process (Jlt)t>o defined by IIj = Yz t for every t > is a continuous Markov chain: 

• its positive transition rates are from a partition tt = to a partition tt® j with 
J C I having at least two elements; the value of such a transition rate is 

1 ( k \ |J| 

k>\J\ { + ' (fc 1 ,...,fc| J |)e(N*)i-'l, ^'-"'W u =i 
ki+...+k\j\=k 

• if tt is a partition of [0, 1] with k non-empty blocks (b±, . . . , bj.) then 

i=l 

Proof. By construction (Y n ) n is a *P([0, l])-valued Markov chain such that for every n £ N, 
Y n +\ = Y n or Y n+ i is a coarser partition of [0, 1] than Y n obtained by merging several blocks 
of Y n in one block. 

Let tt = (bi)i£i be a partition of [0, 1] and let J be a subset of I with at least two elements. Set 
bj =: Uj g j6j. Given the event {Y n = tt}, Y n+ i is equal to tt® j if among the R points U\, . . . , Ur 
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uniformly distributed in the interval [0, 1] at least one point falls in every block 6j for i 6 J 
and none falls in [0, 1] \ bj. Therefore, 

+oo k 

P(Y n+1 = ^ J \Y n = 7r) = J>({*}) E II Leb (^) 

k>\J\ (ii,-,ifc)eW fc (J)tt=l 
+oo 



fc>|J| (fci,-..,fc|j|)e(N*)l J l 

Given the event {Y n = n}, Y n+ i is equal to 7r if Ui, . . . , Ur fall in one block of the partition ir. 
Therefore, 

p(Y n+1 = MYn =«) = j2 k w) ( E Leb (^) fe ) = -} + E mi - • 

fc>2 is/ is/ 

The generator of (IT)t>o is then ^4 = j3 e (Q — I) where Q is the transition matrix of (Y n ) n . 
Let us now compute ¥ no (Ht ^ ir) where ir is a partition of [0,1] with fe non-empty blocks 
(bi, . . . ,bk) coarser than or equal to ttq. {IT h. means that all points that are chosen 
simultaneously in the interval [0, 1] before time t belong to the same block of ir. Therefore, by 
decomposing {IT ^ tt} according to the value of Zt and the number of points falling at the 
same time, we obtain: 

/ +oc (ft t) n n 
K (U t h ir) = exp(-/3 £ t)(l+E E "(W) ■ ■ ■ "({'»}) n(Leb(6i) r '+. . .+Leb(6 fc ) rt ) 

n=l ' ri,...,r n >2 i=l 

We expand the product 

n 

JJ(Leb(6i) r * + ... + Leb(6 fc ) p *)= E Leb(b u jK . .Leb(b Un ) r - 

»=1 (ni,...,«„)G{l,...,fe}" 

and use that 

E -(M) l*< w = i ( - iog ( i - ^> - ^ 

to obtain: 

+0O 



p„(n lt .) = exp(- At )(i + ES E n<-'°s(i-^)-^ 

n=l (ui,...,u„)e{l,. ..,*;}" t=l 

The second sum is also equal to: 

e L.".,Jn(->^-^)-^)"' 

v y i=l 



{m,...,n k )e. 
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This yields 

p„ ( n lt .) _ exp,- Ai) (i + E e nSt-'^-^'-^n 

n=1 (i!i,...,n fc )eN fc ! =1 * 

th\ + . . .-\-rife=7i 

( a+.+VV , n Leb(6i) Leb(6j) \ 
= oqj^- At + ^(-108(1--^)--^)) 

= exp( t log(l-^)- 4 £log(l-^)). 

1=1 

□ 

5.3. Asymptotics for cover time. A vertex x is said to be isolated at time a if no loop 

(n) 

a passes through x. We call 'cover time' the smallest a such that K n has no isolated 
vertex at time a and we denote it T n . 

Let us assume that the intensity of the killing measure is proportional to the size of the graph: 
with e > 0. If we use the algorithm described in Proposition 5.2 to define Ca , then 
we need to choose uniformly at random an average of af3 e E(R) = points in {1, ... , n}. 

If the points are drawn one by one and not by packs of random sizes, then the solution of the 
classical coupon collector's problem provides that the values of T n is around e(l + e)nlog(n) 
for large n. This analogy holds: 

Proposition 5.4. Assume that the intensity of the killing measure on K n is K n = n£ with 
£ > 0. 

For every a £ R, the number of isolated vertices at time a n = e(l + e)n(log(n) + a + o(l)) 
converges in law to the Poisson distribution with parameter e~ a . In particular, W£ ^ e ^ — log(n) 
converges in law to the Gumbel distribution^ . 

Proof. For each vertex x, let I n)X denote the indicator of the event l x is isolated at time a n \ The 
number of isolated vertices at time a n is S n = Ylx=i^nx- For every k £ N*, its k-th factorial 
moment E(S n ) k := E(S n (S n - 1) . . . (S n - k + 1)) is the sum of P{I n ,xi = ... = I n ,x k = 1) 
over all (ordered) /c-tuples of distinct vertices (xi, . . . , x^). As -P(/n,xi = • • • = In,x k = 1) is the 
probability that no loop in VC an intersects the subset F = {x\, . . . , Xk}, 

P{I n ,xi = ■■■ = I n ,x k = 1) = exp ( - a n fj,(£, £ intersects F)^j 

and 

H(£, I intersects F) = fi(VC(X)) - n(VC{F c )) 

= log ( det(G) n x *) - io § ( ^(G {Fc) ) n A -) 

x£{l,...,n} x£F c 

Therefore, E{S n ) k = k\J2 Fe v k ({i,...,n}) ( di^JTLTiT^ ) 

In our setting, for every x € {1, . . . , n}, X x = = n(e + 1) — 1 and for every D C {1, . . . , n}, 
^The cumulative distribution function of the Gumbel distribution is x \-¥ cxp(— exp(— x)). 
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det^))" 1 = (AW + ^PI-^AW + 1 - \D\). Thus 

E(S n ) k = n (n - 1) . . . n - k + 1 (l - n ■ A ( 1 + — 

V n(l + ej / V ne 

= JJ(1 - ±) exp (*(log(n) - ) + fc^L(_ |^ + 4) + 

- LJ - n V nell + e zn z (1 + er e z n^e^ 

i=i 

As a n = ne(l + e)(log(n) + a + o(l)), we deduce that E(S n )k converges to exp(— ka) for every 
k G N*. The convergence to the Poisson distribution with parameter e _a follows from the 
theory of moments. □ 

Remark 5.5. 

(i) The distribution of the number of isolated vertices at a time a can also be obtained by 
the inclusion-exclusion principle; the probability that there exists r isolated vertices at 
time a is 

n ~ r i i i ■ 

Yi-iy (i — r Q(r+j) (i + — -r a - 

> r! j!(n-r- j)! 1 n{l + e)> 1 ne ' 

(ii) In the statement of Proposition 5.4, we assume that the intensity of the killing measure 
K n is proportional to the number of vertices n: K n = ne. In fact the same proof shows 
that Proposition 5.4 also holds if we only assume that \^ n ^ converges to +oo and 
replace e by e n = ^ in the statement. Let us note that if l0 g(* n ) converges to a constant 

c > then E(S n )k converges to nik = exp(— ka + |-) for every k G N* and l —^r- tends 
to +oo for any t > 0. 

5.4. Asymptotics for coalescence time. Let r n denote the coalescence time, i.e. the first 

(n) 

time at which C a has only one block. The following theorem shows that the cover time and 
the coalescent time have the same asymptotic distribution: 

Theorem 5.6. Let us assume that K n = ne with e > 0. Set a n = e(l + e)n(log(n) + a + o(l)) 
where a is a fixed real. 

(n) 

For every k G N, the probability that T>£ a ^ consists only of a component of size greater or 

e ~ak 

equal to 2 and k isolated points converges to exp(— e~ a ) — — as n tends to +oo. 

k\ 

In particular, — — log(n) converges in distribution to the Gumbel distribution. 

ne(l + e) 

Remark 5.7. 

(i) Theorem 5.6 shows that e(e + l)nlog(n) is a sharp threshold function for the connect- 
edness of the random graph process (^i n ^) Q >o defined by the loop sets (D£a ) a . 

(ii) Such properties have been proven for the first time by Erdos and Renyi in [4] for the 
random graph model they have introduced. To facilitate comparison, the following 
theorem states some of their results in a slightly different way from that used in [4] : 



22 



YVES LE JAN AND SOPHIE LEMAIRE 



Theorem (Erdos and Renyi, [4]). Let G(n,N) denote a random graph obtained by 

( N\ 

forming N links between n labelled vertices, each of the , n -. graphs being equally 



\(2)/ 

likely. 

For every c£l and every k £ N, the probability that G(n, |_§(log(n) + c)J) contains a 

connected component of size n — k and k isolated points converges to exp(— e~ c ) ^ as 
n — > tends to +00. 

. „ u(£, £ passes through exactly two vertices) 

111) Let us note that the ratio of j— - , . r-r converges to 

V ' _ n(V£({l,...,n})) 

— ^(e + l + (e + l) 2 log ( f — ) ) 1 as n tends to +00 . The limit is an increasing function 
of e that converges to 1 as e tends to +00. Therefore, it is not surprising to obtain 
properties similar to the Erdos-Renyi model for large e. 

In fact, we can show using the same proof that Theorem 5.6 holds if we replace e by a 
positive sequence (e n ) n such that liminf log(n)e n > 0. 

n— >+oo 

By contrast, if we let e converge to as n tends to +00, the loop sets can have very 
large loops: a result of Y. Chang in [3] implies that if e ~ n~ d with d > 1 then 

n— >+oo 

F(3££VC {n \ covering {l,...,n}) -»• 

log(n) n-> + oo 

The rest of the section is devoted to the proof of Theorem 5.6. 

(n) 

Let A n> k denote the event ir DC a „ consists only of a component of size greater or equal to 2 and 

(n) 

k isolated points', let V n be the number of isolated vertices in VCa^ and let B n be the event 

in) 

l T>Ca„ has at least two components of size greater or equal to 2'. We have 

F(A n>k ) = P(V n = k) — P({V n = k}C\ B n ). 

By Proposition 5.4, F(V n = k) converges to e~ e "^j- • We shall prove that ¥{B n ) converges 
to 0. For a subset F, let qF, n denote the probability that F is a block at time ct n . 

\n/2\ 

r=2 FeV r {{l,...,n}) 

As VCan i s independent of VC^ \ DC 1 ^ , we have 

q F , n = P{DC%f) is connected) P (no loop in VC^ intersects F and F c ). 

For a sufficiently large set F, we shall simply bound qFn from above by 

P(no loop in VC^} intersects F and F c ). 

For small F, the probability F(T>C^ F ^ is connected) is small. Instead of computing it, we 
shall consider its upper bound by the probability that the total length of non-trivial loops 
on VC^a^ are greater or equal to \F\ that is P(X^eF ^""'^ > |^|) where Nx a ' F ^ denotes 
the number of crossings of a vertex x by the set of non-trivial loops included in F at time a. 
Therefore for a small set F, we shall use that 

qF,n < P(^ N^ an ' F) > |F|)P(no loop in VC an intersects F and F c ). 
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We start by stating two lemmas: Lemma 5.8 provides an upper bound for the probability 

(n) 

that no loop in T>CaJ intersects F and F c . Lemma 5.9 gives an exponential inequality for 
P(^] J , ef ivi"™'^ > \F\) based on Markov's inequality. To shorten the notations let V(n) denote 
the power set of {1, . . . , n} and let V r (n) consist of subsets of {1, . . . , n} of cardinality r. 

Lemma 5.8. For every nonempty proper subset F of {1, . . . ,n}, let A n ^p denote the event 'no 
loop in T>Ca] passes through both F and F c ' with a n = e(e + l)n(log(n) + a + o(l)). 
For every S e]0, 1[ there exists ns > such that for every n > n$ and for every n 1 "* 5 < r < ^, 
E F& v r (n) F ( A n,F) < ^eM-^n 1 - 5 log(n)). 

Proof. Let r £ {1, . . . , n/2} and let F be a subset of {1, . . . , n} with r elements. 

/ \ /det(G( F ))det(G( FC ))\«n 
P(A»,f) = exp - a n ^({£, I intersects F and F C })J = ( V det(G) -J . 

In our setting, det(G^) = (n(l + e)y\ D \ +1 (n(l + e) - \D\)^ for every subset D. Therefore, 
F(A n>F ) = (1 + ^^(1 - Using that (») < 1 (^(1 - r)-(-r) ( see for 

example [2], formula 1.5 page 4), we obtain: Y*,Fev r (n) ^( A n,F) < 4^ exp(-n/ n (£)) where f n 
is the function on ]0, 1/2] defined by: 

f n (x) = x log(x) + (1 - x) log(l -x) + u n e(l + e) log(l + ^~^ ) for x e]0, 1/2], 

with u n = log(n) + a + o(l). 

The function f n is of class C 2 in the interval ]0, 1/2] and the first two derivatives of f n are: 



f n (x) = log(x) - log(l -x) + u n 



l-2x 



1+ 



x(l — x) 



• f!l(x) = ® n ^ X ^ xa-x) — where Q n is the polynomial function of second order defined 

by: Q n (y) = i + y(^rb) - (2 + i^K) + y 2 ^) (ipfej + 2u «) for ever y f ■ 

Thus / n has the following properties: 

(i) lim x _, 0+ f n (x) = 0, / n (l/2) > 0, 

(ii) lim x .^ 0+ f'(x) < 0, f(l/2) = 0, 

(iii) f£ = has at most two solutions in ]0, 1/2]. 

A C 2 function in the interval ]0, 1/2] that verifies conditions (i), (ii), (iii) has the following 
property: for every c e]0, 1/2] such that f n (c) is positive the minimum of f n on the interval 
[c, 1/2] lies on c or 7j (see Lemma A. 2 in the appendix). In fact, (f n (l/2)) n converges to +oo, 
whence for sufficiently large n, ini x e[c,i/2] fn(x) = f n (c). 

To conclude, let us consider f n at point —g for some 5 g]0, 1[. 

1. «51og(n) 1 1 / n~ & (l-n- & ) 



fn(zs) = zH- + i 1 ~ Zs) + ^ £ (! + e) log 1 + 



ra° n" n" n" \ e(l + e) 

Therefore f n (^j) is equivalent to log(n). We deduce that for sufficiently large values of n 
and for every x <E [-4-, 1/2], f n (x) > f n (-k) > log(n) which ends the proof. □ 
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Lemma 5.9. Let 5 and 5 be two positive reals such that < 5 < 8 < 1. Assume that K n = ne 
with e > and set a n = e(l + e)n(log(n) + a + o(l)) with a G K /or every n G N. 
There exists n s j > sticfo i/iai /or every n > rigg, and F G 'P(n) 2 < \F\ < n 1 ^ 5 , 



]T > \F\) < 



--I-FI 



Proof. To prove Lemma 5.9, we shall apply the Markov inequality to the random variable 
exp(# ^xeF Nx° n ' FS> ) for a well-chosen positive real 9. 

The generating function of the vector (Nx ,x G F) has been computed for any finite graph 
G = (V,£) in [9] 4 , page 37: 

e ( n s ^ a,F> ) = ( det (h^ + a -( i - s ^ g s] x F )) v (sa;)xev g]0, i]V - 

In our setting, the generating function of J2xeF satisfies: 

E ( s ^FNi a - F> \ = Uet(sI ]Fl + A (n) (l -s)G (F) )Y a VO < s < 1 



where A (n) = n(e+l)-l and G x F y = A( J, +1 (I{i=y} + A (n) + \„| F | ) f° r x,y £ F. The computation 
of that determinant is detailed in the Appendix (Lemma A. 3). We obtain: 

det( S J |F| + AW(1 - .)G^) = (1 - ^ I ) |J?| - 1 (l + (1 " ') AW + r_ 1 | F | )- ( 5 - 4 ) 

Equality (5.4) can be extended to s G M and the determinant is positive for every < s < prpi- 
By the Markov inequality, for every 9 > 0, 

Nt n ' F) > \F\) < exp(-^ F , n (0)) (5.5) 

where ipF,n{6) = — log ^E(exp(— 9\F\ + &^2 xe p N x an ,F ^))\ The expression of the generating 

function of ^2 xe p N^ an ' F ' shows that the value of ipF,n{&) depends on F only via |F|; we denote 
it if}\F\ n (0). The value of Vv,n(#) is 

^ n {9) = 9r + a n ((r - 1) log(l + + log(l - (e s - 1) ^ ~ - _ - ) 

for every > such that e e < ^z~r- 

Set 3 rn = , A< )+1 and # rn , = log(/3 rn ) for every integer 2 < r < n. We shall use inequality 

(5.5) for = 9 r>n . Before, let us note that mm((3 r>n , r G {2, . . . , re 1 " 5 }) converges to +oo and 
max( ^ r, °^ ) ^ , r G {2, . . . , n 1-(5 }) converges to 0. Therefore, for sufficiently large values of n 
and r G {2, . . . ,n 1_<5 }, the two conditions 6> r , n > and exp(0 ri „) < -^-j- are satisfied. 



^Formula (4.3) page 37 in [9] contains a misprint: the exponent 1 in the left-hand side term of the equation 
has to be replaced by a. 
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It remains to bound from below il> r ,n(&r,n) 



AA9r,n) = r(log(AW + 1) - i log(a n ) - ± log(r - 1)) 



+ a n ((r - 1) log(l + ^-4) + log(l - (r - 1) - 1 ) 

V A( n ) + 1 A( n ) + 1 - r 

Using the classical lower bounds 

log(l - t) > -t - t 2 for < t < 1/2 and log(l + 1) > t - -t 2 for t > 

along with max ( (r — 1) f^'" , r G {2, . . . ,n 1 ~ s \ ) — ► 0, we obtain 

V / n— H-oo 

i>r,n(0r,n) > -(log(n) - log(r - 1) - I njl ) - 7„ ;2 ~ I n ,3 

where, for sufficiently large values of n 

• I n>1 = log(l + Tri) + log(log(n) + a + o( l)) < 2 log(log(n)); 

• In,2 = a n (r - l)r (A(n)+ ^-i +1 _ r) < 2r^(log(n) + a + o(l)); 

• /„,3 = a n (r - l)(/3 r , n - 1) 2 ( 2(A( 4 +1) 2 + (A (n) r +i_ r)a ) < 5 + 2(r - 1). 

Therefore there exists a constant M such that for sufficiently large values of n and for every 
2 < r < § , Vv,n(#r,n) > §(log(n) - log(r) - 2 log(log(n)) -M). In particular for 2 < r < n 1 ^ 5 , 
Vv,n(#r,n) > <^§ log(ra) — 2 log(log(n)) — M). This shows that for every < 5 < 5 and sufficiently 
large values of n, min r6 r 2j ... )n i-«i Vv,n(#r,n) > frlog(n) which ends the proof. 
Let us note that a study of Vv,n shows that for every r G {2, . . . , ra 1_<5 }, and sufficiently large 
values of n, ip r)n has a maximum at a point which is equivalent to Q r ^ n as n tends to +oo. □ 

We can now complete the proof of Theorem 5.6. Set S r , n = ^Fe"P r (n) lF,n- We shall prove 
that the upper bound of F(B n ), £i= 2 2J S r , n converges to as n tends to +oo. 
By Lemma 5.8, 

n 

2 I — S 

^ S r ,n < nexp( — n 1_<5 log(n)). 

r=n l — S 

Let us consider the case of subsets F with r < n 1 ^ 5 elements. Using the notations introduced 
in Lemmas 5.8 and 5.9, we have 5V,n < ^ ex P( — n fn(^) — Vv,n(#r,n))- By Lemma 5.9, for every 
5 G]0, <5[, sufficiently large n and r G {2, . . . , n 1-5 }, 

1 / r 5 r \ 

Sy,n < -p exp - n(f n (-) + -- log(n)) . 

V n In J 

Let us study the function / n (x) = / n (a0 + |xlog(n). By computations we obtain that f n (2/n) 
is equivalent to Xlog(n) as n tends to +oo. The study of f n made in the proof of Lemma 
5.8 shows that, for sufficiently large values of n, f n is greater than f n (2/n) in [2/n, 1/2]. Let 
us introduce a real <5 such that < 5 < 6. We have shown that for n large enough and every 
r G {2, . . . ,n 1 ~ 5 }, S r>n < ±rT^ and thus *£r=2 S r,n < n l ~ 5 ~~ s l 2 . 

By taking <5 = 3/4 and 5 = 2/3 for instance we obtain that P{B n ) < ne's™ 174 lo s(") + n~ 1 / 12 
for sufficiently large values of n. This ends the proof of Theorem 5.6. 
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Appendix A. 

Lemma A.l. For a finite set E, let &e denote the set of permutations of E and let <5® E consist 
of permutations of E without fixed point. For two integers r > 2 and 1 < k < r/2, let ^3fc(2, r) 
denote the partitions of {1, ... ,r} with k blocks, each of them having at least two elements. 
For a r x r matrix A and a real a, set Per^(A) = X^o-e6° a ^1^(1)- ■■^r,a(r) where m(o~) 
denotes the number of cycles in a permutation a (Per^A) introduced in [9] page 41, can also 
be defined as the a-permanent of A with the diagonal elements of the matrix set to zero). 
Another expression of PerJ!(A) is 

r/2 k 

Pev° r (A) = Y,a k E II(uh ^ 4(1)^)^)^(3) ■■■4(K 3 |),,(i))- 

k=l 7r=(7ri,...,7r fc )e t Pfc(2,r) j=l 3 ae6 Vj 

Proof. For a finite set B, let G C B be the subset of &b which consists of cycles of length \B\. 
The decomposition of permutations into disjoint cycles entails that 

per? ( A)=xy e n( e n^W) 

k=l ir=(n u ...,n k )£<$ k (2,r)j = l <Tj66^. u67Tj 



Therefore, it remains to prove that for every k > 2 

k 1 

E II i v(«) = jt ^ A u (l),v(2)K{2),u(Z)---A y {k),v{l)- (A.l) 
creS {l,...,fc} 11=1 "6©{l,...,fc} 

Starting from a permutation G ..,&}, we define a cycle of length /c, F{y) such that 
Ili=i A,F(f)(i) = ^i/(i),i/(2)^i/(2),i/(3) • • • A/(fc),i/(i) by setting 



1/(1) ifzv^ 1 (n) = A; 

+ 1) if v~ x {u) <k — l. 



Conversely, starting from a cycle a of length k, we can construct exactly k different permuta- 
tions v\,. . . , v k such that F(vi) = . . . = F(i>k) = a by setting 



i if j = 1 
oJ'-^i) if j€{2 fc} 



This shows equality (A.l), completing the proof. □ 

Lemma A. 2. Let a < b be two reals. Let f be a real function of class C 2 in the interval ]a, b]. 
Assume that: 

• lim x ^ a + f{x) = 0, f(b) > 0, 
. lim^ Q+ f{x) < 0, f{b) < 0, 

• /" = has at most two solutions in ]a, b]. 

Then for every c g]o, b] such that /(c) > 0, the minimum of f on the interval [c,b] lies on c 
or b. 



Proof. The assumptions on / ensure that the only possible configurations are: 
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(i) there exists p &]a, b[ such that / is a decreasing function on ]a, p[ and an increasing 
function on [p, b] ; 

(ii) there exists two reals p_ < p+ in ]a,b[ such that / is a decreasing function on ]a, p-[, 
an increasing function on and a decreasing function on [p + ,b]. 

Since lim a ,_ >a + f(x) = 0, for every x e]a, b[ such that f(x) > 0, we have in both cases 
inf /(«) = min(/(x), /(&)). □ 

Lemma A. 3. Let a and b be two reals and let n be a positive integer. Let I n denote the 
identity matrix of size n and let J n denote the n x n matrix with all entries equal to one. The 
determinant of the matrix M n = bJ n + (a — b)I n is (a — b) n ~ 1 (a + (n — 1)6). 

Proof. First, det(Mi) = a. Let n > 2. If we substract the column n — 1 to the column n and 
expand the determinant along the column n, we obtain: 

det(M n ) = (a — b) det(M n _i) + (a — b) det(i? n _i) where R n denotes the matrix defined by 

\ b if n = 1 

I bJ n + diag(a — b, . . . , a — b, 0) if n > 2. 

By applying the same transformations to R n , we obtain 

det{R n ) = (o - 6) det(i? n _i) = (a - b) n ~ l b. 

The expression of det(M n ) follows by induction on n. □ 
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